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ABSTRACT 

We examine the velocity field of galaxies around the Virgo cluster induced by 
' its overdensity. A sample of 1792 galaxies with distances from the Tip of the 

C"| | Red Giant Branch, the Cepheid luminosity, the SNIa luminosity, the surface 

brightness fluctuation method, and the Tully-Fisher relation has been used to 

o ■ 

■ study the velocity- distance relation in the Virgocentric coordinates. Attention 

was paid to some observational biases affected the Hubble flow around Virgo. 

We estimate the radius of the zero-velocity surface for the Virgo cluster 
to be within (5.0 - 7.5) Mpc corresponding to (17 — 26)° at the mean cluster 



> , 

distance of 17.0 Mpc. In the case of spherical symmetry with cosmological 
■ parameter Q m = 0.24 and the age of the Universe To = l3.7Gyr, it yields the 

| total mass of the Virgo cluster to be within Mt = (2.7 — 8.9) • 10 Mq in 

reasonable agreement with the existing virial mass estimates for the cluster. 

o 
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1 INTRODUCTION 

The gravitational action of mass of a solitary system of galaxies leads to deceleration of the 
local Hubble flow. As a result, the line of average velocity of neighboring galaxies relative to 
the center of a cluster (or a group) deviates from the linear Hubble relation, going to negative 
values at small distances R < R . Here, _R means the radius of the zero-velocity surface, 
which separates the galaxy system against the global cosmic expansion. As it was shown by 
Lynden-Bell (1981) and Sandage (1986), in the simplest case of spherical symmetry with 
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cosmological parameter A = the radius Rq, depends only on the total mass of a group M T 
and the age of the Universe T : 

Mt = (tt 2 /8G) ■ Rq ■ T _2 , (1) 

where G is the gravitational constant. Measuring R via distances and radial velocities of 
galaxies outside the virial radius of the system R V i r , one can determine the total mass of the 
system independent of its virial mass estimate. Note that both the methods deriving mass 
from internal and from external galaxy motions correspond to different linear scales where 
R is roughly 4 times as large as the virial radius. In reality, galaxy groups and clusters do 
not have perfect spherical symmetry, and cosmology with A = is not true. 

Numerous measurements of distances to nearby galaxies obtained recently with the Hub- 
ble Space Telescope (HST) allowed us to investigate the Hubble flow around the Local Group 
and other proximate groups (Karachentsev et al. 2002, Karachentsev & Kashibadze 2006, 
Karachentsev et al. 2006, 2009). The radii Rq obtained from observations for nearby groups 
around: the Milky Way & Andromeda (the Local Group), M81, Centaurus A, Maffei & 
IC 342, NGC 253 (Sculptor filament), and NGC 4736 (Canes Venatici I cloud) are ranged 
within (0.7—1.4 Mpc). The average ratio of total-to- virial masses for these six groups, derived 
from Rq via eq. (1) and from R V i r , turns out to be < M^jM^ >= 0.60 ± 0.15 (Karachent- 
sev, 2005). But as it was noticed by Peirani & Pacheco (2006, 2008) and Karachentsev et al. 
(2007), in a flat universe dominated by dark energy the resulting Mt{Ro) mass is higher than 
that derived from the canonical Lemaitre-Tolman eq. (1). In the "concordant" cosmological 
model with A-term and Q m as a matter component eq. (1) takes a form 

M T =(n 2 /8G)-Rl-H 2 /f 2 (n rn ), (2) 

where 

f(Q m ) = (1 - O^)- 1 - (Q m /2) • (1 - fi m )" 3/2 ■ arccosh[(2/n m ) - 1]. (3) 

Assuming Q m = 0.24 and H = 72 km s -1 Mpc -1 that corresponds to T = 13.7 Gyr (Spergel 
et al. 2007), one can rewrite (2) as 

(M T /M ) O . 24 = 2.12 • 10 12 (R /Mpc) 3 . (4) 
It yields the mass that is 1.5 as large as derived from the classic eq. (1). This correction 
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The observed infall of galaxies towards the Virgo cluster 3 

leads to a good agreement in average between the Rq mass estimates and virial masses for 
the abovementioned galaxy groups. 

For galaxies around the nearest cluster in Virgo, expected velocity deviations from the 
pure Hubble flow (so-called Virgocentric infall) were regarded in dynamical models by Hoff- 
man et al. (1980), Tonry & Davis (1981) and Hoffman & Salpeter (1982). These authors 
note that with the virial mass of the Virgo cluster M v i r ~ 6 • 10 14 Mq, the radius of the 
zero- velocity surface around the cluster amounts to ~ 27°, i.e. the infall zone covers nearly 
1 steradian of the sky. According to Hoffman et al. (1980), the observed decrease of radial 
velocity dispersion within the angular distance = [0 — 24]° from the Virgo center for 228 
galaxies agrees, in the main, with the Virgocentric infall pattern for the cluster mass men- 
tioned above. Tully & Shaya (1984) considered the phenomena of infall of galaxies towards 
Virgo both in the point-mass and distributed-mass models for the cluster with different 
values of the cosmological parameter and the age of the Universe To. Using Tully-Fisher 
distance estimates for 19 galaxies inside the virial radius of ~ 6° and 14 galaxies outside it, 
the authors ascertain the expected infall with Rq ~ 28°. 

Later, Tonry et al. (2000, 2001) developed a model of the Virgocentric flow basing on 
accurate distance measurements for 300 E and SO galaxies from surface brightness fluctua- 
tions. Their model fits well the observational data on galaxy distances and radial velocities 
for the Virgo cluster distance of 17.0 Mpc and its total mass of 7 ■ 1O 14 M . According to 
their model, our Local Group has a peculiar velocity of 139 km s -1 directed towards the 
Virgo center. 

Teerikorpi et al. (1992), Ekholm et al. (1999) and Ekholm et al. (2000) examined the 
Virgocentric flow with different models of density distribution in the cluster and infered 
expected relations between velocities and distances of galaxies relative to the cluster center. 
Using Cepheid distances to 23 galaxies and Tully-Fisher distances to 96 galaxies, the authors 
conclude that the radius of the zero-velocity surface ranges from 20° to 31°, and the total 
cluster mass is equal to (1 -v- 2) of its virial value. 

During the last decade, the observational database on distances to galaxies in a wide 
vicinity of the Virgo cluster has grown significantly, allowing to determine Rq and, therefore, 
the total mass of the Virgo cluster with better accuracy. 
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2 THE OBSERVATIONAL SAMPLES 

To examine phenomena of the Virgocentric flow, we used distance moduli of galaxies from 
different publications preferring more precise measurements. The main data sources are 
listed below. 

A) Taking luminosity of the Tip of the Red Giant Branch (TRGB) as a standard is the 
most efficient and the most universal method to determine distances to nearby galaxies, as it 
is practically independent of their morphological type. Being applied to galaxy images in two 
ore more photometric bands obtained with WFPC2 or ACS cameras at the HST, the TRGB 
method yields an accuracy of distance measurements of ~ 7%, as it was founded by Rizzi 
et al. (2007). A consolidated list of distances for the Local Volume galaxies is presented in 
the Catalog of Neighboring Galaxies (= CNG, Karachentsev et al. 2004). The CNG sample 
of 451 galaxies has been collected based on two conditions: galaxy distance D < 10 Mpc, if 
a galaxy has individual distance estimate; otherwise galaxy radial velocity with respect to 
the Local Group Vlg < 550 km s -1 . Below we use from CNG the galaxies with only TRGB 
or Cepheid distances, supplying them with new TRGB distances from recent publications 
(Karachentsev et al. 2006, Tully et al. 2006). 

B) The surface brightness fluctuation method, SBF, applying to early type galaxies, 
assumes that the old stellar population (RGB) is prevailing in a total luminosity, and the 
galaxy structure does not suffer from irregularities due to dust clouds. Using this approach, 
Tonry et al. (2001) determined SBF distances to 300 E and SO galaxies with a typical errors 
of ~ 12%. This sample is distributed over the whole sky extending to cz ~ 4000 km s- 1 with 
a median velocity of 1800 km s _1 . 

C) Mei et al.(2007) undertook a two-color ACS/HST imaging survey for 100 early-type 
galaxies situated in the Virgo cluster core (the ACSVCS project). They derived precise SBF 
distances to 84 E, SO galaxies with a typical error of 8%, and revealed a 3D-shape of the 
Virgo cluster to be a slightly triaxial ellipsoid with axis ratios of (1:0.7:0.5). We expanded 
the ACSVCS sample with other precise SBF and TRGB distance measurements in the Virgo 
core made by Neilsen & Tsvetanov (2000) and Caldwell (2006) with ACS/HST and by Jerjen 
et al. (2004) with the VLT. That yields us the total "ACSVCS+" sample of 116 galaxies. 

D) In wide vicinity of the Virgo cluster there are 22 galaxies with distances measured by 
Tonry et al. (2003) via SNIa. This sample is small but has the distance error of only 5%. 

E) Based on the "2MASS Selected Flat Galaxy Catalog" (Mitronova et al. 2004 = 
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2MFGC), Kashibadze (2008) determined distances to 402 spiral edge-on galaxies with radial 
velocities < 3000 km s -1 . A multiparametric NIR Tully-Fisher relation was applied to them, 
yielding a typical distance error of ~ 20%. The zero point of the luminosity - line width 
relation was calibrated by 15 galaxies with Cepheid and TRGB distance measures. 

F) Finally, the former samples of galaxies were supplemented with a compilation of 
distances by Tully et al. (2008, 2009) that have been obtained from optical (B, R or I 
band) one-parametric Tully-Fisher relations. This compilation relays on numerous Hl-line 
and photometric observations carried out by Methewson Sz Ford (1996), Haynes et al.(1999), 
Tully & Pierce (2000), Koribalski et al. (2004), Springob et al. (2005), Theureau et al. (2006) 
and other authors. Zero points of the data were recalibrated by a set of 40 galaxies with 
known Cepheid and TRGB distances. As a last step, we used also distances from a very 
big and important SFI++ sample (Springob et al. 2007), which have not been included in 
Tully et al. (2009) compilation. At total, we used distance estimates for 941 spiral galaxies 
whose radial velocities were limited by 3000 km s . The typical distance error for them is 
~ 20%, though there are some cases with much higher errors due to uncertainties of galaxy 
inclination, presence of interacting companions or HI profiles of low quality. 

A substantial overlap between the two last luminosity-linewidth samples provides con- 
firmation that their zero points are the same and gives rms agreement per measure of 0.40 
mag. As it is seen from Fig. 1 by Tully et al. (2008), there is an excellent agreement in dis- 
tance moduli between the luminosity-linewidth and other (Cepheid, TRGB, SBF and SNIa) 
measures. In particular, for 12 galaxies in our list with both TF and TRGB or SNIa moduli 
the mean distance difference is (—1.4 ± 1.2) Mpc, while for 20 galaxies with SBF or TRGB 
moduli the average difference is only (0.2 ± 0.2) Mpc. 

Following the previous authors (Tully & Shaya 1984, Ekholm et al. 2000) we have formed 
a composite sample of galaxies limiting their angular separation from the Virgo center to 
< 30°. We have considered the radiogalaxy Virgo A = NGC 4486 to be the physical center 
of the cluster as its position is close to the center of X-ray emitting gas. The total number 
of galaxies in this cone volume with apex angle G < 30° is 630 (the 2D sample). 

Limiting the angular separation of galaxies from the Virgo center introduces some selec- 
tion effects into the Virgocentric flow analysis. For this reason we have used also another 
way to form the observational sample, considering galaxies with spatial distances from the 
Virgo cluster center R vc < 30 Mpc (the 3D sample). This approach suffers a drawback too 
because galaxy distances are measured with errors and their significance is different at the 
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proximate and the distant boundary of the spherical volume (the so-called Malmquist bias). 
The total number of galaxies in our 3D sample amounts to 1792, and the fractions of diverse 
subsamples differ significantly from that in the 2D sample. 

Table 1 presents the summary of observational data that we have used. The first column 
indicates subsample kinds, the second one gives typical distance errors expressed in magni- 
tudes. Columns (3) and (5) contain numbers of galaxies in the cone (2D) or in the spherical 
(3D) volume. The sample goodness G, defined as G = (iV/lOO) 1 / 2 ■ a - " 1 , is a useful parameter 
which characterizes a statistical weight of a certain sample (Kudrya et al. 2003). Goodness 
values are indicated in columns (4) and (6). For example, the subsample of galaxies with 
SNIa distances is scanty but its statistical significance is comparable with that of other 
samples because of higher accuracy of distance measurements. As one can see, the galaxy 
subsample ACSVCS+ has the maximum statistical weight in the 2D set, however almost all 
these galaxies are concentrated within the virial radius. In the 3D set the highest goodness 
corresponds to the TRGB sample, but its majority is crowded on the nearby side of the 
examined volume. The last two TF samples exhibit a significant increase in number going 
from the 2D to the 3D samples that is caused by the well-known effect of morphological 
segregation of late-type vs. early-type galaxies along the cluster radius. 



3 RADIUS OF THE ZERO- VELOCITY SURFACE R 

The virial radius of the Virgo cluster R V i r = 1.8 Mpc (Hoffman et al. 1980) corresponds 
to its angular scale of 6.0°, assuming the average distance to the cluster members to be 
17.0 Mpc. Radial velocities and distances relative to the Local Group (LG) centroid for 
259 galaxies in this zone are represented in the top panel of Fig. 1. Here, precise distances 
for most the galaxies were obtained within the special survey ACSVCS with HST (Mei et 
al. 2007). The Virgo cluster members, located in the distance range from 14 to 20 Mpc, 
demonstrate a radial velocity scatter from —800 up to +2300 km s _1 . Foreground galaxies 
are scarcely presented on the panel while background objects tend to lie below the linear 
Hubble regression with the global Hubble parameter H = 72 km s _1 Mpc -1 (Spergel et al. 
2007), showing thereby the expected effect of infall into the Virgo cluster from the opposite 
side. The centroid of galaxies forming the "virial column" at [17.0 ± 1.8] Mpc, marked by 
vertical lines, has a mean velocity +1004±70 km s _1 versus the expected Hubble velocity 
of +1224 km s _1 at the distance of 17.0 Mpc, which can be explained by a peculiar motion 
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The observed infall of galaxies towards the Virgo cluster 7 

of the Local Group ~ 220 ± 70 km s _1 directed towards Virgo. The dotted and solid S- 
shaped curves correspond to a Hubble flow perturbed by a point-like mass of 2.7 ■ 10 14 Mq 
and 8.9 • 10 14 Mq (as the limiting cases discussed below) for the line-of-sight passing exactly 
through the cluster center. 

The distributions of radial velocities and distances for remaining galaxies of the 2D 
sample in close surroundings of Virgo, (6° < < 15°) and in a distant periphery (15° < 
< 30°), are shown in the middle and the bottom panels of Fig. 1. Here, the solid and 
dotted S-shaped lines having lower amplitudes describe the behavior of perturbed Hubble 
flow at angular distances equal 6° and 15°, respectively. These panels display some signs 
of the infall effect too, however, in front of the Virgo, the expected infall is seen barely. 

Considering a set of such kind Hubble diagrams with their different amplitudes of S- 
shaped waves decreasing with the angular distance 0, one can find the quantity of the 
cluster mass that fits the observed infall pattern in the best way. But this approach seems 
us to be not transperent enough. In order to determine Ro and the total cluster mass via it, 
we have converted our observational data into distances and velocities expressed relative to 
the cluster center. The top panel of Fig. 2 shows the layout of a galaxy (G) relative to the 
observer (LG) and the cluster center (C) with angular separation from the cluster center. 
The spatial distance of the galaxy from the center therefore is 

R 2 vc = R 2 g + R 2 c - 2R g R c cos 0. (5) 

Assuming that the galaxy and the cluster center are involved in an almost unperturbed 
Hubble flow (Hf case) with negligible peculiar velocities, we can state the mutual velocity 
difference between G and C in projection onto the stright line connecting them as 

V vc = V g cos A — V c cos /I, (6) 

where \i = A + and tan A = R c smQ/(R g — R c cos 0). 

The distribution of galaxies in the Virgocentric reference frame {V vc , R vc } is represented 
in Fig. 3. The only 391 galaxies with angles A obeying A < 45° or A > 135° from the whole 
2D sample are shown here. Selecting galaxies situated approximately in front and behind 
the cluster is meant to reduce the role of tangential velocity componens which remain still 
unknown. The polygon curve traces the running median with a window of 1 Mpc. The 
median follows roughly the linear Hubble regression with H = 72 km s -1 Mpc -1 (the 
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inclined dashed line) at middle Virgocentric distances of ~15 Mpc, but tends to deviate 
from the H ■ R vc line at smaller scales, crossing the zero-velocity line at Rq ~ 6 Mpc. The 
behaviour of the running median at large scales is strongly skewed by a selection effect due 
to the adopted limit for galaxy velocities Vlg < 3000 km s _1 . 

Another approach can be used also for converting the observational radial velocity of a 
galaxy V g into its Virgocentric velocity (the case of pure Virgocentric flow, Vf). If the galaxy 
is not involved in the general cosmological expansion but is falling instead towards the Virgo 
cluster with a velocity V in ( Fig. 2b), than its radial velocity relative to the observer will be 
expressed as 

V g = V c ■ cos 9 — V in cos A, (7) 
and the infall velocity Vi n itself can be written as 

V in = (V c cos0- V s )secA, (8) 

where the angles G and A are shown in the Fig. 2b. Evidently, the discrepancy between these 
two extreme approaches, Hf and Vf, decreases when A tends to or to 180°. 

As can be easily seen, selecting galaxies both in the cone with apex angle = 30° and 
the angle A entails a loss of galaxies at large Virgocentric distances which becomes apparent 
in the top right corner of the Fig. 3. A selection of galaxies in the volume R vc < 30 Mpc (3D 
sample) reduces the bias appreciably. Fig. 4 represents the sky distribution of galaxies with 
known Virgocentric distances up to 30 Mpc in equatorial coordinates. The galaxies of this 
3D sample are marked as circles and their diameters indicate three distance ranges: (0-12) 
Mpc, (12-22) Mpc, and more than 22 Mpc from the observer. This map exhibits that the 
selection of galaxies by their angle < 30° brings some systematic skews dependent on 
distance. In particular foreground Virgo galaxies are preferentially removed by this angular 
selection. 

The Hubble diagram for 1792 galaxies of the 3D sample is represented in Fig. 5. The 
symbols for objects from different sources of distance data are the same as for Fig. 1. The 
number of galaxies in the 3D sample is roughly 3 times as large as in the 2D sample. It is 
noteworthy, they populate the crucial regions in front and behind the Virgo cluster more 
thoroughly giving us an opportunity for more detailed analyses of the Virgocentric infall. 
The velocity - distance relation for these galaxies with respect to the cluster center is shown 
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in Fig. 6. Its top panel corresponds to the assumption of pure Hubble flow of galaxies (Hf 
case) while the bottom panel represents the case of pure radial motions towards the Virgo 
center (Vf). As previously, the galaxies located far away from the line of sight crossing the 
cluster center, i.e. with 45° < A < 135°, are eliminated in order to reduce the role of unknown 
tangential velocity components. This condition diminishes the number of sampled galaxies 
by 42%. 

Comparison the bottom panel of Fig. 6 with the top one shows that switching from the Hf 
case to the Vf case does not lead to any dramatic changes in the Hubble flow pattern given 
in the Virgocentric coordinates. Some galaxies move along the vertical axis appreciably but 
the total behaviour of the running medians traces the infall of galaxies towards the cluster 
in a similar way. The asymptotic tendency of the median at large distances R vc looks much 
more regular for the 3D sample than for the 2D one. 

Our elimination of galaxies with 45° < A < 135° is a little arbitrary. To estimate the 
response of the Virgocentric flow pattern to changing this condition we have imposed also a 
more rigid constraint, eliminating galaxies with 30° < A < 150°. The corresponding diagrams 
for Hf and Vf cases are represented in two panels of Fig. 7. As it is seen, the more severe 
selection of galaxies by A reduces their number by more than a half. However, the behaviour 
of the running median is almost the same. 

In the bottom panel of Fig. 7 every galaxy from the samples A,B,C, and D (see Table 
1) having a distance error within 12% is supplied by error bars indicating where the galaxy 
should be situated if its distance from the observer changes by ±1(7d. As expected, these error 
bars are much longer for galaxies situated behind the cluster. Changing a galaxy distance 
by ±12% leads in some regions of the {V vc , R vc } diagram to significant displacement of the 
galaxy along the both Virgocentric coordinates, and therefore to appreciable galaxy skips 
relative to the zero- velocity line. (We do not include the longer error bars for the Tully-Fisher 
distances). 

To quantify the uncertainties on the running median curves plotted in Figs. 6 and 7, we 
generated extensive sets of bootstrap realisations. Their results permit us to estimate the 
radius of the zero velocity surface and also its rms error presented in Table 2. The table 
columns contain: (1) size of the window for a running median, taken to be 0.8, 1.0, and 1.2 
Mpc; (2,3) the mean radius of the zero velocity surface Ro and its rms scatter obtained with 
regarding to the assumptions on the pure Hubble flow (Hf) and the pure Virgocentric flow 
(Vf), respectively; here only four samples: A,B,C, and D with precise distance moduli were 
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taken into account; (4,5) the same quantities for the samples E and F with Tully-Fisher 
distances; (6,7) the radii Rq and their errors for the total set of available data on galaxy 
distances. We shell discuss the interpretation of these different estimates of Rq in the next 
section. 

4 DISCUSSION AND CONCLUSIONS 

We have analysed available observational data on distances and radial velocities of galaxies 
in wide surroundings of the Virgo cluster in order to study the Virgocentric infall. The 
main purpose of the present paper is to determine the radius of the zero- velocity surface Rq 
that separates the cluster against the global cosmic expansion. By using this observational 
quantity we are able to estimate the total mass of the Virgo concentrated within R and 
compare it with the virial mass estimates corresponding to roughly 4 times lower scale, R V i r . 
Based on the Hubble diagrams transformed into Virgocentric coordinates and the results 
of our bootstrap numerical experiments given in Table 2, we derive i?o to be in the range 
of 5.0 — 7.5 Mpc with a typical random error of 1.0 Mpc. The derived value of Rq can be 
affected by some systematic circumstances. 

Analysis of the Hubble diagrams in the Virgocentric coordinates suffers the lack of data 
on tangential velocity components of the galaxies. We tried to overcome this drawback with 
assumptions regarding a dominant type of galaxy motions in the Virgo proximity. Conversion 
of the observed radial velocities of galaxies into their Virgocentric velocities was performed 
by us under two extreme kinematic assumptions: almost unperturbed Hubble flow (Hf) or 
almost pure radial flow towards Virgo {Vf). We found that adopting one or another scheme 
doesn't change significantly a general pattern of the Virgocentric infall. Calculating velocities 
in the VJcase yields, at the average, some smaller values of V vc , which causes a slightly larger 
(+0.7 Mpc) value of .Ro- 

The Rq quantities presented in Table 2 tend to be correlated with the smoothing window 
size. When the window changes from 0.8 Mpc to 1.2 Mpc, the radius dicreases on 0.3 Mpc 
at average. Variation the window size in wider range: from 0.5 Mpc to 2.0 Mpc still leaves 
i?o within its random errors. 

The most noticeable systematic variations in the radius Rq are seen in dependence on 
galaxy samples. The use of samples of galaxies with only precise distance measurements (via 
Cepheids, TRGB, SBF, and SNIa) yields the radius R within 6.5 — 7.5 Mpc, while the use of 
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less precise Tully-Fisher distances leads to Ro estimates around 5.2 Mpc. The physical origin 
of this difference is clear. Probably, the instance is just related with the fact that a typical 
distance error for Tully-Fisher method (~ 20%) corresponds at the Virgo distance to a linear 
scale of 3.4 Mpc that comparable with the virial diameter of the cluster (3.6 Mpc). Large 
random errors can throw galaxies over the "virial column" diluting the shape of S-wave infall. 
Nevertheless, we should stress that two completely independent sets of observational data 
on distances to early-type and late-type galaxies lead to quite compatible values of R . As 
can be seen from Table 1, two samples with Tully-Fisher distances (E+F) are about 2 times 
as large in galaxy number as the samples (A+B+C+D) with precise distances. Because we 
estimate Ro based on the running median without regard to galaxy weights (distance errors), 
the radius Ro for the total sample turns out to be closer to that for the former samples. 

Assuming Ro for the Virgo cluster in the range of (5.0 — 7.5) Mpc, we obtain from 
eq. (4) the total mass of the cluster to be within = (2.7 — 8.9) ■ 10 14 Mq. This value is 
consistent with the virial mass estimates for the Virgo: (5.5 ± 1.4) • 10 14 Mq (Hoffman et 
al. 1980), (7.5 ± 1.5) ■ 1O 14 M (Tully & Shaya 1984) and 7.0 ■ 1O 14 M (Tonry et al. 2000) 
normalized to the same Hubble parameter Ho = 72 km s _1 Mpc -1 . It should be reminded 
that this agreement is on very different scales as the zero velocity radius, Ro, is roughly 4 
times as large as the virial radius. This means our results are consistent with there being 
no additional mass outside the virial radius of the Virgo cluster. 

Two items are worth mentioning in conclusions. The stated method of estimating Rq 
value from observational data relays on the assumptions that a cluster has spherically sym- 
metric shape and that galaxy motions around the cluster are regular without a significant 
tangential component. Deviations from the simple spherical infall scenario have been consid- 
ered both theoretically and observationally by many authors (Peebles 1980, Davis & Peebles 
1983, Lilje et al. 1986, Tonry et al. 2000, Mould et al. 2000, Watkins et al. 2009). In partic- 
ular, Tonry et al. (2000) studied a peculiar velocity field around the Local Group and Virgo 
cluster in the presence of second attractor (Great Attractor) situated in Hydra-Centaurus. 
As seen from their Figs. 20,21, the Great Attractor with a total mass of 9 • 10 15 Mq situated 
at a distance of 43 Mpc generates significant distortions in the velocity field over the area 
of ~1 steradian. Besides, Tully (1988) found so called effect of "Local Velocity Anomaly" 
arising by a push of the Local Void. According to Tully et al. (2009), a bulk motion of ~260 
km s -1 toward the Supergalactic pole may be attributed to the Local Void. Recently, cosmic 
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velocity flows in the Local Universe were also studied by Erdogdu et al. (2006), Haugbolle 
et al. (2007) and Lavaux (2009). 

We would emphasize that observational abilities to determine Rq and Mt with a better 
accuracy have not been exhausted yet. With the derived Rq around 7 Mpc, the zero-velocity 
surface is situated at a distance of about 10 Mpc from the observer, i.e. on the Local 
volume edge. Hence, one can expect to find some number of the LV galaxies with radial 
velocities ~ (1000 — 2000) km s^ 1 , residing in front of Virgo. A new updated version of CNG 
(Karachentsev et al. 2010) contains about 40 candidate dwarf galaxies, like VCC1675, with 
appropriate radial velocities and rough distance estimates. Their precise distances could be 
easily measured via TRGB with facilities of ACS at HST in the one-object- per-one-orbit 
mode. The addition of more galaxies with measured distances in the range 9 — 11 Mpc from 
the Local Group would allow a better measurement of Rq, in particular these galaxies could 
be used to better rule out the case where there is a large amount of mass outside the virial 
radius. 
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Table 1. Distance datasets for galaxies in/around Virgo. 



Sample a m N(2D) G(2D) N(3D) G(3D) 



A 


TRGB + Ceph 


0.15 


37 


4.1 


264 


10.8 


B 


SBF (Tonry) 


0.25 


81 


3.6 


189 


5.5 


C 


ACSVCS+ 


0.17 


113 


6.2 


113 


6.3 


D 


SNIa (Tonry) 


0.10 


13 


3.6 


22 


4.7 


E 


TF (IR) 


0.40 


56 


1.9 


260 


4.0 


F 


TF (opt) 


0.40 


330 


4.5 


941 


7.7 



Table 2. Radius Ro and its rms error (in Mpc) obtained from bootstrap realisations. 



Window Primary distances TF distances All 



(Mpc) (samples A+B+C+D) (samples E+F) samples 





Hf 


Vf 


Hf 


Vf 


Hf 


Vf 


0.8 


6.79 


7.78 


5.07 


5.58 


5.15 


5.82 




±1.12 


±1.13 


±1.04 


±1.13 


±0.96 


±1.04 


1.0 


6.59 


7.54 


4.87 


5.38 


5.00 


5.65 




±1.05 


±1.03 


±0.88 


±0.97 


±1.00 


±1.07 


1.2 


6.46 


7.42 


4.71 


5.24 


4.80 


5.52 




±1.01 


±0.97 


±0.78 


±0.87 


±0.95 


±1.03 
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Figure 1. The radial velocity vs. distance relation for galaxies in the Virgo cluster region with respect to the Local Group 
centroid. Galaxy samples with distances derived by different methods are marked by different symbols. The inclined line traces 
Hubble relation with the global Hubble parameter Hq = 72 km s — 1 Mpc -1 . The vertical dashed lines outline the virial zone. 
Two S-shaped lines correspond to a Hubble flow perturbed by virial masses of 2.7 ■ 1O 14 M (dotted) and 8.9 • 10 14 M Q (solid) 
as the limiting cases within confidence range. Top: the cluster core within angular distance < 6°, both the S-shaped lines 
indicate the expected infall at = 0°. Middle: wider surroundings with 6° < < 15°, the S-shaped lines indicate the infall 
at = 6°. The typical distance error bars for each dataset are shown. Bottom: peripheric regions with 15° < < 30°, the 
S-shaped lines indicate the infall at = 15°. 
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Figure 3. Hubble flow around Virgo in the Virgocentric coordinates for 391 galaxies with < 30° and small angles A regarding 
to the line of sight. The inclined dashed line corresponds to Ho = 72 km s — 1 Mpc -1 , and the polygon curve traces the running 
median with a window of 1 Mpc. The dotted and solid curves correspond to a Hubble flow perturbed by a point-like mass of 
2.7 ■ 1O 14 M0 and 8.9 ■ 1O 14 M0 respectively as the limiting cases. 
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Figure 4. Sky distribution of 1792 galaxies with Virgocentric distances R vc < 30 Mpc in equatorial coordinates. Circles of 
different diameters indicate three distance ranges from the LG. 



© 2010 RAS, MNRAS 000, []]-?? 



20 I. D. Karachentsev, O. G. Nasonova (Kashibadze) 



+ x X X + x _ 

x +A ++ , X + ++ ^ 

+ A X + t X *W* « + + 

+ + x A++ x 

^.x + xj W *x 

+v - X A „ X. X + 1 



X»x + 



,4 



ti t ,>^ + >J % 4, . ^ jfr.^ + . -K + ^fcX.^' f + t, 



» + A * 



4* 




A Af A x 



A 

w 



i 

10 



15 

R LG , Mpc 



_J 

20 



• TRGB + Ceph 

& SBF (Tonry) 

▼ ACSVCS+ 

■ SNIa (Tonry) 

x TF (IR) 

+ TF (opt) 

3D sample (N = 1792 gg) 

I 

25 



30 



Figure 5. The Hubble diagram for 1792 galaxies with Virgocentric distances less than 30 Mpc. Symbol indications the same 
as in Fig. 1. 
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Figure 6. The Hubble flow in the Virgocentric coordinates for 1045 galaxies with R vc < 30 Mpc. Top: the case of almost pure 

Hubble flow. Bottom: the case of almost pure Virgocentric infall. Only the galaxies with A < 45° or A > 135° arc presented. 

The dotted and solid curvesjsprrespond to a Hubble flow perturbed by a point-like mass of 2.7 ■ 10 Mr?, and 8.9 ■ 10 14 M,-r, 
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Figure 7. The same Hubble diagram as in Fig. 6, but only for 467 galaxies with the angles A < 30° or A > 150°. In the bottom 

panel all galaxies with precise distances are supplied with error bars indicating where the galaxy should be situated if their 

distance from the observer changes by ±lcrn. _ 
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